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Abstract. We express AKNS hierarchies, admitting reductions to matrix NLS and matrix 
mKdV hierarchies, in terms of a bidifferential graded algebra. Application of a universal 
result in this framework quickly generates an infinite family of exact solutions, including e.g. 
the matrix solitons in the focusing NLS case. Exploiting a general Miura transformation, we 
^ , recover the generalized Heisenberg magnet hierarchy and establish a corresponding solution 

formula for it. Simply by exchanging the roles of the two derivations of the bidifferential 
graded algebra, we recover "negative flows" , leading to an extension of the respective hierar- 
chy. In this way we also meet a matrix and vector version of the short pulse equation and 
also the sine-Gordon equation. For these equations corresponding solution formulas are also 
derived. In all these cases the solutions are parametrized in terms of matrix data that have 
to satisfy a certain Sylvester equation. 
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1 Introduction 

•l-H 

X 

$_i ■ A unification of some integrability aspects and solution generating techniques has recently been 

achieved for a wide class of "integrable" partial differential or difference equations (PDEs) in 
the framework of bidifferential graded algebras pQ . The hurdle to take is to find a bidifferential 
calculus (i.e. bidifferential graded algebra) associated with the respective PDE. In particular, 
a surprisingly simple result (Theorem 3.1 in pQ and Theorem[T]below) then generates a (typically 
large) class of exact solutions. This has been elaborated in detail for matrix NLS systems in 
a recent work [2]. The present work extends some of these results to a corresponding hierarchy 
and moreover to related hierarchies. It demonstrates how to deal with whole hierarchies instead 
of only single equations or systems in the bidifferential calculus approach and shows moreover 
that certain relations between hierarchies find a nice explanation in this framework. Except for 
certain specializations, we deal with "non-commutative equations", i.e. we treat the dependent 
variables as non-commutative matrices, and the solution formulas that we present respect this 
fact. 

In Section [2] we introduce some basic structures needed in the sequel. Section [3] presents 
a bidifferential calculus for a matrix AKNS hierarchy. We derive a class of solutions of the 
latter and address some reductions of the hierarchy. In Sections H] and [5] we show that, in the 
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bidifferential calculus framework, a "reciprocal" [3] or "negative" extension of the hierarchy 
naturally appears. "Negative flows" have been considered previously via negative powers of 
a recursion operator (see e.g. [U El [7] and also [8j [91 [TOj [11] for other aspects). In our 
picture, these rather emerge as "mixed equations" , bridging between the ordinary hierarchy and 
a "purely negative" counterpart. 

Section [6] elaborates this program for a "dual hierarchy". Here we recover in the bidifferential 
calculus framework in particular a well-known duality or gauge equivalence between the (matrix) 
NLS and (generalized) Heisenberg magnet hierarchies [I2l [T3], [141 [I5l [THJ [TTl UHl HH] - Section [7] 
contains some concluding remarks. 



2 Basic structures 

Definition 1. A graded algebra is an associative algebra O over C with a direct sum decompo- 
sition 



n = o r 



r>0 

into a subalgebra A = fi and ^4-bimodules f2 r , such that 

n r n* c n r+s . 

Definition 2. A bidifferential calculus (or bidifferential graded algebra) is a graded algebra O 
equipped with two (C-linear) graded derivations d,d : U — > Q, of degree one (hence dU r C Q r+ , 
dffl C r2 r+1 ), with the properties 

dod = 0, dod = 0, dod + dod = 0, 
and the graded Leibniz rule 

d(xx') = (dx)x' + (-l)W, d(xx') = (dx)x' + (-l) r xdY, 
for all x S O r and x' £ ^> 

For any algebra A, a corresponding graded algebra is given by 

n = A®/\(c N ), (2.i) 

where f\(C N ) denotes the exterior algebra of C , N > 1. Defining graded derivations d, d on A, 
they extend in an obvious way to Q such that the Leibniz rule holds and elements of f\(C N ) are 
treated as constants with respect to d and d. 

Given a bidifferential calculus, it turns out that the equation 

dd<f> = d<f)A d(f>, where 6 A, (2.2) 

has various integrability properties [TJ [2]. By choosing a suitable bidifferential calculus, this 
equation covers in particular the familiar selfdual Yang-Mills equation (in one of its gauge- 
reduced potential versions), but also e.g. discrete integrable equations [D[2]. In the next section 
we demonstrate that, by choosing an appropriate bidifferential calculus, (|2,2p reproduces matrix 
AKNS hierarchies. 

The (modified) Miura transformation 

[dg - (dg)A]g- 1 = dcP, (2.3) 
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where dA = (dA)A, is a hetero-Backlund transformation between (|2.2p and the dual equatior^ 

d{[dg-(dg)A]g~ 1 ) = 0. (2.4) 

In the present work we concentrate on the case where A = 0. We note that (j2.3|) and (|2.4p are 
equivalent if d has trivial cohomology. But this is in general not the case. 

Exchanging d and d in (|2.2p . we get a different equation. Dealing with hierarchies, such an 
exchange leads to what we call the reciprocal hierarchy. If A = 0, exchanging d and d in (|2.4p . 
simply amounts to replacing g by g~ l . 

3 AKNS hierarchies 

Let Bq be the algebra of complex smooth functions of independent variables t±, t2, £3, • • -, B an 
extension by certain operators (specified below), and A = M&t(m,m,B), m > 1. Let m = 
m i + m 2 with m,j G N, and let V be the projection 

T> — T> — I Omi XJT12 1 (n -1 \ 

^ - Hmi,m 2 ) "In ' ^ A,1 > 

\ K, rri2 X m 1 "tn2Xm2 / 

where / m denotes the m x m identity matrix. If the dimension is obvious from the context, we 
will simply denote it by /. It will also be convenient to introduce the matrix 

3.1 NLS system 

A particular bidifferential calculus on A is determined by 

df = [V, f]Ci + [Vd x J]C 2 , df = f x d + \[d t2 + dl /]C 2 , (3.2) 

where Ci, C2 is a basis of /\ 1 (C 2 ), and we set x = t\. Here S is the extension of Bq by the partial 
derivative operator d x . Evaluation of (|2.2p yields 

[P,0t a ] = {P,0 a « B } + 2(^ a! )[P,0]-2[P,0]^ > (3.3) 
using the familiar notation for commutator and anti-commutator. The block-decomposition 

<P = j( P q ) = ( P q ), (3-4) 
\ q p J \ -q -p J 

(where p, p,Q, q have size mi x mi, mi x 7712, m\ x 777-2, ^2 x wi, respectively), results in the 



NLS system^ 

gt 2 = qxx - 2qqq, qt 2 = -q xx + 2qqq, (3.5) 

together with 

Px = -qq, (3.6) 

where we set a "constant" of integration to zero. As a consequence of the form of V, there is no 
equation for p. Though at this point we could simply set it to zero, this would be inconsistent 
with further methods used in this work (cf. Remark [3|). 



introducing B — [dg — (dg)A]g~ 1 , this equation reads dB — 0, and as a consequence, taking also dA = (dA)A 
into account, we find that dB — B A B. If, as in familiar cases, this (partial) zero curvature equation implies 
B = (dg^g'" 1 , then (|2.4[) is gauge-equivalent to d[(d<7')<;' _1 ] = 0, so that the term involving A in (|2.4[) can be 
generated by a gauge transformation. It is nevertheless helpful to consider the modified equation (|2.4|) in order 
to accommodate more easily certain examples of integrable equations in this formalism [T|. 

2 This matrix NLS system apparently first appeared in [20]. See also the list of references in [2] and in 
addition [21] [22] [23 [24] . 
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3.2 Extension to a hierarchy 

Another bidifferential calculus on A is determined by 

df=lVE x J}( 1 + [VE ll J}( 2 , d/ = A~ 1 [lEA,/]Ci + ^ 1 [IE^/]C2. (3.7) 

Here Ea and E^ are commuting invertible operators, which also commute with V, and B is the 
extension of Bq by these operators. Introducing /m = E^/E^" 1 and f~\\] = E^" 1 /E^, (|2.2|) results 
in 

(A-i/ - V cj> + 4>_ [x] v) _ M - ^ + 

= - + _ [A] (A- 1 / - 7> + 0_ [A] P) . (3.8) 

In the following, Ea will be chosen as the Miwa shift operator, hence 

/ ±[A] (*!, i 2 , t 3 , ...) = f(h ±\,t 2 ± A 2 /2, h ± A 3 /3, • • •), 

with an arbitrary constant A (see e.g. [25] ) . The generating equatiorH (|3.8p already appeared 
in [25] and we recall some consequences from this reference. Expanding (j3.8|) in powers of 
the arbitrary constants (indeterminates) A and [J,, we recover (j3.3|) as the coefficient of A 1 // . 
Decomposing the matrix <j) into blocks according to (|3.4p . this results in 



(A 1 l-p + p _ [x] )_ M (fj, i/_ p + p_ M ) 

= (/^~ 1 -^-p + p-m)„ [a] (A _1 i"-p + p_[A]) + (qq)-[\], (3.9) 

and 

A _1 (q - g_ [A] ) +p_[ A ]9 = ^"H 1 ? - (3-10) 

A _1 (<Z[a] - <z) + gp[A] = - <z) + (3-11) 

Again, there is no equation for p. Since the last equations separate with respect to A and 
they imply 

A -1 (<? - g_[A]) -Qx-{p- P-[\])q = 0, 

A _1 (Q-9-[A]) - ?-[A]^b + ?-[A](P-P-[A]) = °- ( 3 - 12 ) 

Multiplying the first equation from the right by Q-[\], the second from the left by q, using (|3.6|) 
and adding the resulting equations, we obtain 

(P-P-[X] +Mq-[\])x = \p-P-[\] + Mq-[\],P ~ P-[\]]- (3-13) 
If we set integration constants to zero, this leads to 

x- 1 (p-p-[\]) = -qq-[\]- (3-14) 

Using this equation, (|3.12|) becomes 

A _1 ((? - q-[\]) ~ q x + Mq~[X]Q = 0, A _1 (g [A ] - q) - q x - \qqyx\q = 0, (3.15) 



3 In some publications such an equation has been called a "functional representation" of the corresponding 
hierarchy, see [25] and the references cited therein. It seems to be more appropriate to call it a "generating 
(differential) equation" . 
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which are generating equations for a hierarchy that contains the (matrix) NLS system. Together 
with (|3.6p . this leads to 

(A- 1 /-p + p_ [A] )_ [M] (/i _1 /-p + p_ M ) - (pJ-M 

= i^ 1 ! ~ P + P-[ll\)-[\\{^~ 1 I -P + P-{\\) ~ (Px)-[\], (3.16) 



which is a generating equation for the potential KP hierarchy |26t I25) 4 I . If we think of p as 
determined via (I3.6D in terms of q and q, then the last equation is a consequence of (|3.15p . 
Furthermore, the two equations (|3. 12[) are the linear respectively adjoint linear system of the 
KP hierarchy (cf. [27]) in the form of generating equations. 
Let us recall that 



E\ = exp 



E I = E A " S ^) where B = { d ^\d t2 , l -d t „ . . .V (3.17) 

v n>l J n>0 ^ ' 



and s n are the elementary Schur polynomials. Expanding (|3.15|) in powers of A thus leads to 
s n (-d)(q) - qs n _ 2 (-d){q)q = 0, 

s n (B)(q) - qs n _ 2 (B)(q)q = 0, n = 2,3,.... (3.18) 



For n = 2 we recover f|3. 5 [) . For n = 3, and after elimination of ^-derivatives using (|3.5p . we 
obtain the system 

qt 3 = qxxx - 3(q x qq + qqq x ), qt 3 = q xxx - 3(q x qq + qqq x ), (3.19) 

which admits reductions to matrix KdV and matrix mKdV equations (see also e.g. [2lJ[22j[2E])- 
In the same way, any pair of equations in (|3.18p can be expressed in the form 

qt n = Qn(q,q,q x ,q x , ■ ■ ■ , q x ™,q x n ), qt n = Qn(q, q, q x , q x , ■ ■ ■ , q x n ,q x n ), (3.20) 

by use of the equations for q tk , q tk , with k = 2, . . . ,n — 1. For n = 4, we find 

qti = qxxxx ~ ^qqq xx - ^q xx qq - 2qq xx q - Gq x qq x - 2qq x q x - 2q x q x q + Gqqqqq, 

qti = ~q xxxx 

+ 4qqq xx + Aq xx qq + 2qq xx q + 6q x qq x + 2qq x q x + 2q x q x q - Qqqqqq. 
Remark 1. Expansion of (|3.7p leads to 

df=J2 A m [Ps m (5),/]Ci +^2» n [Ps n (d)J}( 2 , 

m>0 n>0 
m>l n>l 

Selecting the terms with the same powers of A and fi, this suggests to define 

d( m ,„)/ = [Ps m {d)J}Ci + [Vs n (d),f]( 2 , d( m ,n)f = [sm+i(<9),/]Ci + [s n+ i(d),f]( 2 . 

For any choice of non-negative integers m, n, this determines a bidifferential calculus. With 
m = and n = 1, we recover (|3.2p . 



4 Its first member is the potential KP equation (4p t3 — p xxx — 6(p x ) ) x — Pt 2 t 2 + 6{Px,Pt 2 ] = 0. 
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Remark 2. (|2.2p is the integr ability condition of the linear equation 

dlr = (d<f>)V + v&m (3.21) 

for an m x m matrix ^, where v is a constant (cf. pQ). If ^ is invertible, this is (|2 . 3[) with 
A = z^J. Evaluation for the above bidifferential calculus leads to 

- tf_ [A] ) = (7>0 - 0_ [A] P)¥ + v(V^> - *_ [A] P). 
In terms of ^ given by 

this takes the form 

- </>-[A]) = ("P + Tty - 

which is a generating equation for all Lax pairs of the hierarchy. Expanding in powers of A, the 
first two members of this family of linear equations are 

ip x = Lip, ip t2 = Mip, 

where 



L = vV + [V,< 



v q 
q 

M = v 2 V + y[P, 4>] + [V,4>] 2 + {V,4 )x }=( ^ + m t 2Px m _ + qx 

\ vq — q x qq 

This constitutes a Lax pair for the NLS system (|3.5p . In order to obtain a more common Lax 
pair for the NLS system, we have to eliminate p x via (I3.6P and add a constant times the identity 
matrix to L and M, together with a redefinition of the "spectral parameter". See also [2]. 



3.3 A class of solutions 

So far we defined a bidifferential calculus on Mat(m, m, B). In the following we need to extend 
it to a larger algebra. The space of all matrices over B, with size greater or equal to that of 
no x hq matrices, 

Mat no (£)= Mat(n',n,B), 

n',n>no 

attains the structure of a complex algebra A with the usual matrix product extended trivially by 
setting AB = whenever the sizes of A and B do not match. For the example introduced in the 
preceding subsection, we set no = 2. Let VL be the corresponding graded algebra (|2. 1 j) . For each 
n > 2 and a split n = n\ + 112, we choose a projection matrix Vr niin2 \ of the form (13. ip . Then we 
can extend the bidifferential calculus defined in (j3.7|) to A by simply defining the commutators 
appearing there appropriately, e.g. for an n x m matrix / we set 

[VE X , f] = V (nun2) E x f - fV {mum2) E x . 

For the following general result, see [H[2]. 
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Theorem 1. Let (fi,d,d) be a bidifferential calculus with Q = A® f\(C N ) and A = Mat no (B), 
for some uq G N. For fixed n > uq, let X, Y G Mat(n, n, B) be solutions of the linear equations 

dX = (dX)S, dY = (dY)S, 

and 

RX -XS = -QY, Q = VU, (3.22) 

with d- and d-constant matrices S,R G Mat(n, n, B), U G Mat(m, n, B), V G M.&t{n,m,B). 
If X is invertible, then 

$ = UYX~ l V G Mat(m, m, B) (3.23) 

satisfies 

d<t) = (d0)(f> + d0, (3.24) 
with some m x m matrix 9. By application of d, this then implies that (ft solves (|2.2p . 



Now we apply this theorem to the bidifferential calculus associated with the hierarchy intro- 
duced in Section |51 We fix m, 112, and write y for Vr nijV , 2 )- The linear equation dX = (dX)S 
is then equivalent to 

\-\x - X_ [A] ) = (VX - X_ [X] V)S. 

d and d-constancy of S means that S is constant in the usual sense (i.e. does not depend on the 
independent variables ii,i2>---) an d satisfies ["P,£] = 0, which restricts it to a block-diagonal 
matrix, i.e. S = block-diag(S'i, S2). Decomposing X = Xj + X Q into a block-diagonal and an 
off-block-diagonal part, and using [P,X(i\ = 0, [P,X Q ] = JX D = —X Q J with 

we obtain 

(X d - X d> _ [x] )(I - XV S) = 0, X (I - XVS) = X 0| _ [A] (I - XV S), 

where V denotes the projection V — J ■ The first equation implies that X d is constant. We 
write Xd = A d - Noting that VS and VS commute, the solution of the second equation is 

X = A e E fc>i ( ^ s)fc * fc e- E <>i (ps) '\ 

with a constant off-block-diagonal matrix A Q , hence 

X = A d + A Q S where 3 = e"^ J , {(S) = ^S*t t . 

k>l 

A corresponding expression holds for Y, 

Y = B d + B Q S. 
Now (|3.22p splits into the two parts 

RA d + VUB d = A d S, RA + VUB = A D S. 
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Assuming that A d is invertible, we can solve the first of these equations for R and use the 
resulting formula to eliminate R from the second. This results in 

SK-KS = VU, (3.25) 

where 

K = -A d 1 A , U = U{B + B d K), V = AfV. (3.26) 
Note that K = K Q and JU = —UJ (whereas JU = UJ). Next we evaluate (|3.23D . 

4> = UYX- l V = U(B d + B S)(A d + AoSy'V = (UB d + UB 3){I - KB^V 
= {UB d + (U - UB d K)a){I - KS^V = US(I - KS)- l V + UB d V. 

Using the identity (J — K3)" 1 = ( I + KS)(I — (KS) 2 )" 1 , this decomposes into 

<Pd = <Pd,o + USKS(I - (KSfyW, d , o = UB d V, 

and 

<P = UE(I-(KS) 2 )- 1 V. 
All this leads to the following result, which generalizes Proposition 5.1 in [2]. 
Proposition 1. Let 



(3.27) 



(3.28) 





S 


s 


U 


U 


V 


V 


size 


m x m 


n 2 x n 2 


m\ x ri2 


ni2 x m 


ni x mi 


ni x ni2 



be constant complex matrices, and let K {of size n\ x 112) and K [of size ri2 x n\) be solutions 
of the Sylvester equations 



SK + KS = VU, 



SK + KS = VU. 



Then 



q = UE(I n2 - KZKZy'V, q = UE(I m - KEK1 



(3.29) 



(3.30) 



where 



-as) 



t(s) = ^2s k t k , 



k>l 



solve the hierarchy (|3.15p . Furthermore, 

p = UEKE(I ni - KEKE^V 
solves (|3.14p and then also the potential KP hierarchy ()3.16|) . 



(3.31) 



Proof. The expressions (13^9]) . (I3T301) and (I3T3TT) follow, respectively, from (l3T25l . (13381) and 
(|3.27p . by writing 



K 



K 
K 



S 0_ 
-S 



u 



o_ u 
-u 



V 



V 
V 



, (3.32) 
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S = block-diag(H, H), and using (|3,4p . From Theorem Q] we know that (|3,30p and (|3.3ip solve 
the hierarchy equations (|3.9|) . (|3.1U|) and (|3.11|) . It should be noticed, however, that on the way 
to the hierarchy (|3.15p the step to (|3.14[) involved a restriction. Hence we have to verify that p 
given by (|3.3ip actually solves (|3.14p . Noting that 

H_ [A] =Z(I-\S)-\ E_ [X] =E(I + XS), ZEK = KEZ, 

where I stands for the respective identity matrix, Z = H _1 — KEK and Z = H _1 — KEK, we 
have 

Z_ [x] = E-\l - XS) - KE(I + XS)K = Z- X(E~ 1 S + KESK) 

= Z - X(E~ l S + KE[VU - KS}) = Z(I - XS - XZ^KEVU) 
= Z(I -XS- XEKZ~ l VU), 

where we used the second of equations (I3.29p . Now we find that p = UEKZ~ 1 V satisfies 

p-p-ix] = ue(kz~ 1 z ax] - {i + xs)k)zzI x] v 

= Ue(R(I - XS - XEKZ^VU) - (/ + XS)K\ZZ\ X ^ 



XUE [KEKZ~ l VU + SK + Ks)zz 1 



A] 



V 



--VU 



= -XUE(KEK + Z)Z- 1 VUZZ\ X] V = -Agg_ [A] . ■ 

There are matrix data for which the Sylvester equations (|3.29p have no solution. But if S 
and —S have no common eigenvalue, they admit a solution, irrespective of the right hand side, 
and this solution is then unique (see e.g. Theorem 4.4.6 in [29]). 

Remark 3. Evaluated for the bidifferential calculus (13. 7|) . (13.24j) reads 

x-\<j> - 4>- { x]) = (Vd> - <f>-[x]P)<f> + re- e_ [x] v. 

Using (|3.4p and a corresponding block-decomposition for 6, 

s r 
-f —s 

this equation splits into the system 



A" 


-\p- 


P- 


-w) = 


(P 


- P-[\])P ~ m + * - «-[A] 


A" 


-\<1- 


Q- 


-[A]) = 


(P 


-P-[\])q- QP + r, 


A" 


-\<i- 


Q- 


-[A]) = 


-Q 


r _ [A ]P-f_[A], 


A" 


-\p- 


V- 


-[A]) = 


~(i 


'-[X]Q, 



which implies q x = —qp + r, q x = —qp — f, and p x = —qq. With the help of these equations 
we recover (|3.12p and thus also (|3.13p . We observe that the solution generating method based 
on Theorem [J also imposes differential equations on p. The corresponding generating equation 
is actually the direct counterpart of (|3.14p . which the solutions determined by Proposition [1] 
satisfy. 
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Remark 4. If 9 in (|3.24p is not restricted and if we have trivial d-cohomology, then (|3.24p 
is equivalent to (|2.2|) . But d given by (|3.7|) has non-trivial cohomology. Indeed, the following 
1-form is d-closed but not d-exact, 

,=«»(■<? :)<^{y °)c, 

Here a, b only depend on the parameter A, respectively fi, and a star stands for an arbitrary entry 
(not restricted in the dependence on the independent variables and the respective parameter). 
Adding p on the r.h.s. of (I3.24p would achieve equivalence with (12. 2|) . 

3.4 Reductions 

Let us consider the substitution dt 2k —dt 2k , k = 1, 2, in the hierarchy equations (|3. 18[) . 
By inspection of ()3.17j) . it implies s n (d) h-» (— l) n s n {— d), hence maps f)3. 18[) into 

s n (d)(q) - qs n - 2 (d)(q)q = 0, 

s n (-d)(q) - qs n _ 2 (-d){q)q = 0, re = 2, 3, ... , 

which has the same effect as 

q h-» eg", q h-> eg w , 

where e = ±1, and w is the involution on the algebra of matrices either given by the identity 
map or by complex conjugation (u) = *), or the anti-involution either given by transposition 
(uj = j) or by Hermitian conjugation (uj = f) (see also [2] ji- 
lt follows that the odd-time part of the hierarchy, expressed in the form (|3.20p . is consistent 
with the reduction condition 

q = e<f, (3.33) 

which reduces any of its pairs to a single member. In particular, (|3.19p becomes the matrix 
mKdV equation 

qt 3 = qxxx - ^{q x q u q + qq u qx), e = ±1. 

The reduced hierarchy is therefore a matrix mKdV hierarchy. 
After the replacement 

t 2 k ^ i*2fe, A; = 1,2,..., 

with i = v— T, so that dt 2k i— )■ —idt 2k , also the even-time equations of the hierarchy are consistent 
with the above reduction (|3.33[) . provided we choose for uj complex conjugation (*) or Hermitian 
conjugation (f), so that i w = — i. Then (|3.5p becomes the matrix NLS equation 

iq t2 = -q xx + 2eqq u) q, e = ±1. 

This matrix version of the NLS equation apparently first appeared in [20] , The corresponding 
reduced hierarchy is a matrix NLS hierarchy. 

Since Proposition Q] provides a class of solutions of the original hierarchy in terms of matrix 
data, we should address the question what kind of constraints a reduction imposes on the latter. 

5 Note that this may require restricting the size of the matrices. For example, if lu is complex conjugation, we 
need mi = m%. 
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3.4.1 Reduction using an involution 

If u> is one of the involutions specified above, setting 

5 = 5^, U = ee'U UJ , V = e'V LU , K = eK u , (3.34) 

with e' = dbl, and arranging that 5 = H w , i.e. £(— 5 W ) = £(— S) = — £(»S) W , achieves the reduction 
condition (|3.33p . This forces us to set 

mi = 7712, n\ = n 2 . 

Renaming m; to m and rij to n, this leads to the following consequence of Proposition [TJ 

Proposition 2. Let S,U,V be constant n x n, m x n, respectively n x m matrices, and let K 
be a solution of the Sylvester equation 

SK + KS" = VU. 

Then 

q = ±UE UJ (I n — eK u EKE u )~ 1 V u , where E = e~ 5(5) 

with 

' ^S 2fc+1 i 2 fe+i if w = id ("real" mKdV), 



as) 



k>0 



^5 2fc+1 t 2fc+ i + i^5 2fc t 2fc if uj = * (NLS-mKdV), 



k>0 k>l 

solves the m x m matrix "real" mKdV, respectively NLS-mKdV hierarchy. 

If the involution is given by complex conjugation, in the focusing NLS case the solutions 
obtained from this proposition include matrix (multiple) solitons. See |2] for corresponding 
results for the respective NLS equation, the first member of the hierarchjlj. 

3.4.2 Reduction using an anti-involution 

In this case the reduction condition (I3.33P can be implemented on the solutions determined by 
Proposition [1] by setting 

S = S U , U = V U ', U = eV", K w = K, K" = K, 

and arranging again that E = E u . For the anti-involutions specified above, we are forced to set 
n i = which we rename to n. Then we have the following result. 

Proposition 3. Let S, V , V be constant matrices of size nxn, nx mi and n x ni2, respectively. 
Let K, K be (with respect to uj) Hermitian solutions of the Sylvester equation^ 

SK + KS U = VV 10 , S"K + KS = VV U . 

Then 

q = V u E u (L n - eKEKE u )- l V, where E = e~^ s) 



In [2] we used a bidifferential calculus different from the one chosen in the present work. But the resulting 
expressions for exact solutions are the same. 

7 Here the previous matrix K has been redefined with a factor e. 
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as) 



with 

' ^S 2fc+1 t 2 fc+i if oo = T ("real" mKdV), 

^5 2fc+1 t 2fc+1 + i^5 2fc t 2fc if w = t {NLS-mKdV}, 

fc>0 fc>l 

solves the mi x m 2 matrix "real" mKdV, respectively NLS-mKdV hierarchy. 

If the involution is given by Hermitian conjugation, in the focusing NLS case the solutions 
obtained from the last proposition include matrix (multiple) solitons. Corresponding results for 
the respective NLS equation, the first member of the hierarchy, have been obtained in [2]. 

4 The reciprocal AKNS hierarchy 

Exchanging the roles of d and d in (13.Tf) . we have 

d/ = A" 1 [E A ,/]Ci + ^ 1 [^,/]C 2 , d/ = [PE A ,/]Ci + [PE M ,/]C 2 . (4.1) 

Here the Miwa shift operator is defined in terms of a new set of independent variables, ii, 
i = 1,2,..., and / £ Mat(m, m, B), where B is the algebra of smooth functions of these variables, 
extended by the Miwa shifts. Now (j2.2j) results in 



= (p-^^-^-p^^P-A- 1 ^-^)). (4.2) 
In terms of 

ip = 4> — xV, 
where x = t±, the above equation reduces to 

(<P ~ fa - <P-fr]) = (<P- ¥>_ [At] )_[A](¥> - <P-[\])- (4.3) 
Expanding this in powers of A and /i, to order A 2 ^ it yields 

(<Px 2 )x-[<Ps,<Pi a ] = °- ( 44 ) 

This is nothing but the nonlinear part of the potential KP equation. More generally, (I4.3D is the 

nonlinear part of the potential KP hierarchy as obtained from (|3.16j 1. To order A 3 /i, (|4. 3j) yield 



s 



X 



Remark 5. (|4.4p is related to the generalized Heisenberg magnet (gHM) equation 25f 2 = 
[S,Sxx] (see e.g. [13] [30], EU [H]) for an m x m matrix S as follows. The gHM equation is 
compatible with the constraint 5 2 = /. Writing S = 2(p x — I, the constraint reads ip x 2 = ip x , 
which implies 'Px'-Pxxfx = 0, and the gHM equation becomes <*pi 2 = [<fxifxx\-> setting a constant 
of integration to zero. Then [y>s,y>j 2 ] = (<Ar 2 )z> hence tp satisfies the nonlinear part of the 
potential KP equation. 



The nonlinear part can be extracted via the scaling limit e — > of the potential KP hierarchy, assuming 
p(ti,t2, . . .) = ip(ii,t2, . . .) for the KP variable p, where i„ = e n t n . 
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Remark 6. The general linear equation ()3.21| ^l leads to 



Setting 

\ n>l 

this takes the form 

which is a generating equation for all Lax pairs of the reciprocal hierarchy. The first two 
members of this family of linear equations are 



= Lip, 



ih. = Mip, 



where 



L 



v 1 ip s = v 



1 - Px -Qx 

q x Px 



M = v~ 2 ^ 2 - i/-Vt 2 



4 + 2p 2 . + p s 2 - q s q s ) - v 1 pt 2 v 2 {q s + p x Qx ~ QxPx) ~ v 1 



-v 2 {q x + qsPx ~ PxQx) + v x qi 



v 2 (Pi 



q x qx) + v l pt 



4.1 A class of solutions 

We apply again Theorem [TJ Using (|4.ip . dX = {dX)S takes the form 

x-'(x - x_n) = (rx - x_ { - x] r)s-\ 

assuming that S is invertible. Decomposition as in Section 13.31 leads to 



X = A d + A Q S, Y = B d + B D S 



where 3 



k>i 



These are the same formulas we obtained in Section l3.3[ but with S replaced by S -1 . From 
(|3.22p we obtain, however, the same Sylvester equation, SK — KS = VU, using the same 
definitions as in (I3.26p . Furthermore, we obtain again (I3.27P and (j3.28|) . and thus the following 
counterpart of Proposition [TJ 

Proposition 4. Let 





S 


S 


U 


U 


V 


V 


size 


m x m 


n 2 x n 2 


mi x n 2 


m 2 x m 


ni x mi 


n 2 x m 2 



be constant complex matrices, where S and S are invertible, and let K {of size ni x n 2 ) and K 
{of size n 2 x m) be solutions of the Sylvester equations 



SK + KS = VU, 



SK + KS = VU. 



9 Using the definitions of Section T3, 2 1 this linear equation reads d 1 !/ = (d(f>)'i> + v&fy , since in the present section 
we exchanged the definitions of d and d relative to those of Section 13.21 
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Then ip = (f>— xV ', where (f) is given by (j3.4j) with the components 

q = UE(I n2 - KEKE)' 1 ?, q = UE(I ni - KEKE)~ l V, 

p = UEKE(I ni - KEKE^V, p = UEKE{I n2 - KEKE)~ X V , 

and 

E = e~^\ S = «S) = £5"***, 

fc>i 

solves the hierarchy (|4,3ft . 

5 The combined hierarchy 

The bidifferential calculus determined by 

df = [7>E Al ,/]Ci + [VE X2 J]C 2 + + / ^2 1 [E At2 ,/]C 2 , 

df = X^[E Xl ,f]Ci + A2 1 [E A2 ,/]C 2 + [VE^JKi + [PE M2 ,/]C 2 , (5.1) 

contains (|3.7p and ()4. 1|) . and thus combines the corresponding hierarchies. Here £i> (2, Ci> C2 
is a basis of /\ 1 (C 4 ), Aj and fii are indeterminates, and / £ Mat(m, m, ,6), where i3 is now the 
algebra of complex smooth functions of independent variables t±, t 2 , £3, . . ., ti,i 2 ,---, extended 
by the Miwa shift operators. Then ([22]) is equivalent to (pTKI) . (IQ) . anoS 

(V - iT\<t> ~ 0-fl)-[A]) (A^ 1 / - - 0-[A]^)) 

To order A , respectively this yields 

A- 1 ^ - 0_ [A] ) S = (P - - 0_ [A] P) - (^0 - 0-[A]^)(^ " (5-2) 

M -1 ^-^])* = (^-^- 1 (0-^ ] ))[P,0] - [P,0]„ [ - ] (P-^ 1 (0-0_ [ - ] )). (5.3) 

To order A°/x° we have 

4>xx = [r -^[r, 

Using (|3.4p . this becomes 

fe + Px? + qpx = q, q xx + P x q + qps = q, (5.4) 

and p xS = -(qq)x, Pxs = ~(qq)x, which integrates to 

Px = ~qq, Px = ~qq, (5.5) 

setting constants of integration to zero. 
To order A, (\5.2h leads to 

¥t2X = ~ [fx, {P, <Px} + [P, if} 2 } , 

10 This equation alone is obtained from (|2.2[) using the bidifferential calculus determined by df — ["PEa,/]C + 
/i _1 [E M , /]£ and df — A _1 [Ea, f]C + fPE^, /]£. The remaining equations of the hierarchy can be recovered from 
the linear system with this (restricted) bidifferential calculus. 
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which decomposes into 

Pt 2 x = - [ps , qq] - q x q~x + q x q x , 

Pt 2 x = \px, qq] + q x q x - q x q x , 

qt 2 x = (ps - l)q x + qxPs + q x qq + qqq x , 

qt 2 x = -q x (px - i) - q x p x - q x qq - qqqx- 

Furthermore, to order fi (|5.3|) yields 

ft 2 x = -[fhiVM] - [v,<px]}, 

which leads to 

qi 2 x = -Pt 2 q - qpt 2 - Psqs + qsPx, qi 2X = -Pt 2 q - m 2 - Psqx + qsp x - 

In the same way as in the preceding sections, we arrive at the following result. 
Proposition 5. Let 





S 


s 


U 


U 


V 


V 


size 


m x m 


n 2 x n 2 


mi x ri2 


rri2 x n% 


ni x mi 


n 2 x m 2 



be constant complex matrices, where S and S are invertible, and let K {of size n\ x n 2 ) and K 
{of size n-2 x n{) be solutions of the Sylvester equations 

SK + KS = VU, SK + KS = VU. 

Then <f> given by fj3.4|) with the components 

q = UE{I n2 - KEKS)- 1 ?, q = UE{I ni - KEKE)~ 1 V, 

p = UEKE{I ni - KEK2)- X V, p = UEKE{I n , 2 - KEKE)^, 

where 

5 = e-t( s \ E = e^V, f(5) = ^ 8% + ^ S~%, 

k>l k>l 

solves all equations of the combined hierarchy. 

This simply extends Propositions [1] and U] by adding the respective expressions for £(S). 

5.1 A reduction 

Let q, q, p, p be square matrices of the same size. Setting 

q = eq, p = p, where e = ±1, (5-6) 
the system (|5.4p . (|5.5p reduces to 

Px = ~eq 2 , q xx = -(I - 2p s )q + -q(I - 2p^j. 

This reduction corresponds to the choice w = id in Section 13.41 and extends more generally to 
the "odd flows" of the combined hierarchy. Imposing the conditions (|3.34|) on the matrix data 
of Proposition [51 then leads to the following result. 
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Proposition 6. Let S, U , V be constant n x n, m x n, respectively n x m matrices, and let K 
be a solution of the Sylvester equation SK + KS = VU . Then 

q = ±UE(l n - e(iifE) 2 )~V, p = eUEKE(l n - e(KE) 2 yW, 

where 

fc>i fc>i 
sofee ifoe odd pari of the combined hierarchy with the reduction condition (15.61) . 

5.1.1 Short pulse equation 

Let us impose the additional condition that p is a scalar times the identity matrix. Then we 
have 

Px = -eg 2 , Qxx = (I ~ 2px)q- 
Writing 

P = -Ax-z)I m , (5.7) 

with a new dependent scalar variable z, the last system is turned into 

z x I = 2eq 2 , q xx = z x q. (5.8) 
In terms of u(x, z) given by 

u(x, z(x, x)) = 2q(x, x), (5-9) 
we obtain 



2Qx =u x + z x u z = u x + -u u z , z x u = 2q xx = z x [u x + -u u x 

The change of independent variables requires z x ^ 0. The last equation is then equivalent to 

Uxz = u - ^(u 2 u z ) z , (5.10) 

which is a matrix version of the short pulse equation. The latter apparently first appeared in [32] 
(see also |33} 134] ) and was later derived as an approximation for the propagation of ultra-short 
pulses in nonlinear media [35] . It was further studied in particular in [361 EH ESI EHl HOI HTj H21 
143] , Of course, we have to take the additional condition into account that u 2 has to be a scalar 
times the identity matrix. This is achieved if 



U = 

i=l 

where e^ej + e^ej = 2r]ijl with rjij = iSij (Clifford algebra), since then u 2 = (u,u)I, where 
u = (ui, . . . , u r ) J and (u, u) = Yli=i r Hj u i u j- I n this case f|5. lOj) becomes 

u xz = u-U(u,u)u z ) z . (5.11) 

This vector version of the short pulse equation is different from those considered in \42\ [4"4"j . In 
the following example we obtain an infinite set of exact solutions of the 2-component system 
via Proposition [6l 
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Figure 1. Parametric plots (with parameter x) of the two components u\ and U2 for a 2-soliton solution 
of the 2-component short pulse equation (|5.11|) with e = — 1, determined by the data of Example Q] with 
N = 2, si = 1, s 2 = 2 and C = diag(l + 3i, 1 - 3i, 2 + i/2, 2 - i/2) at x = -3, -1, 1. 



Example 1. We can alternatively express the solution given in Proposition [6] in the a priori 
more redundant form 



q = ±UE(l n - e(KE) 2 ) V, p = eUEKE(l n - e(KE) 2 ) V, 



2\-h 



with 



SK + KS = VU, 



CI 



where C is any constant n x n matrix that commutes with S. More precisely, the introduction 
of C allows us to fix some of the freedom in the choice of the coefficients of the matrices U, V. 
The following choices involve further restrictions, however. We consider the case m = 2, n = 2N, 
and choose 



S = diag(si 7 2 , • • • , s N I 2 ), U = ( o\ ... a\ ) , V 



( h 



where o\ is the respective Pauli matrix. The Sylvester equation is then solved by 
1 



A" 



ij 



Si + s 



■o-i, 



i,j = 1,...,N. 



Choosing C block-diagonal where the 2x2 blocks on the diagonal are a linear combination 
of I2 and the Pauli matrix 03 = diag(l, —1), then C commutes with S. Furthermore, it follows 
that KE consists of 2 x 2 blocks (KE)ij which are linear combinations of the off-diagonal Pauli 
matrices o\ and 02- It further follows that ((KE) 2 )jj is diagonal. As a consequence, the inverse 
of I2N — e(AH) 2 also consists of diagonal 2x2 blocka 11 !. Because U consists of off-diagonal 2x2 
blocks, we conclude that q given by the above formula is an off-diagonal 2x2 matrix. Hence 
its square is proportional to the identity matrix I 2 . Since N G N is arbitrary, we thus have an 
infinite family of exact solutions of the system (|5.1ip with r = 2, where the components of the 
vector u = (ui,U2) J are given by u = u\a\ + U2&2- Fig. [T] shows a plot of a 2-soliton solution. 

Remark 7. In order to obtain a Lax pair for the short pulse equation, we start with 



fa 



L —I ) ib, 



11 This is quite obvious if we regard the matrices as N x iV matrices over the commutative algebra of diagonal 
2x2 matrices. 
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with L and L taken from Remarks [2] and [6l respectively. Without imposing a reduction, the 
integrability conditions are (|5.4|) and (|5.5|) (modulo an integration with respect to x). Writing 

i/j(x, x) = x( x -> z ( x i x))-, we fhid 



it v r\ l u u \ e 2 

Xx = { L--I) X -z xXz = -[ u _\ x --u Xz 



I 


-u z 


eu z 


-I 


\{ 


V 




eu 



¥~t u2 \u + ^u 2 u z 



where we applied the reduction conditions (15.6P and used (15. 7]) , (15.8|) and (15. 9|) . Using the 
symmetry it i— > —u of the short pulse equation, we recover the Lax pair given in [36] for the 
scalar case and with e = — 1. 



6 Dual AKNS hierarchies 

For the bidifferential calculus determined by (|3.7p . the Miura transformation (j2.3[) with A = 
takes the form 

^{g-g-ixfig- 1 =V<f>-</>- [x] V, (6.1) 

and the dual equation (|2.4p becomes 

M _1 - (ff-M'T 1 )-^) = A" 1 (^-[a]^" 1 - (^-[a]^ _1 )-m^) • 

Multiplying from the left by g~^_r , and from the right by g, this can be written as 

Z^ 1 (sl^ - GCf A] P*)- M ) = A ~* (^m^ " (^M^)-[A]) ■ (6-2) 
To order fjP this yields 

A" 1 {g~ l Vg - (g~ l Vg) AX] ) = (gZj x] Vg) x . (6.3) 
Applying a Miwa shift with —[/j] and subtracting the result from this equation, we obtain 

^ (gZ\ x {Pg-{9Z\ x] VgU^ x 

= >r x H- 1 (g- x Vg-{g- x 'P9)-[x\ ~ QT^-M + (^ _1 ^)-[a]-m)- 

Since the r.h.s. is symmetric in A, /x, this implies (|6.2p up to an ^-independent "constant of 
integration". Hence (|6.2p reduces to (|6.3p . The first non-trivial equation resulting from an 
expansion of (|6.3p in powers of the indeterminate A is obtained as the term linear in A, 

(g-'Vgh, - (g- l Vg) xx = -2 ((g- 1 ) x Vg) x . (6.4) 

6.1 Miura transformation and relation 

with the generalized Heisenberg magnet model 

Introducing 



S = g-'Jg 
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(see also [IS]), we have the identities 

S tn = [S,g- 1 g tn ], n = l,2,..., (6.5) 

and 

S 2 = I => S X S = -SS X =► S X 2 = -\{S,S XX } => (SS X ) X = 1 -[S,S XX ]. (6.6) 



The dual hierarchy equation (16. 4p can now be expressed as follows, 

S t2 =S xx + 2(g- 1 g x (I + S)) x . 

In order to express this solely in terms of S we use the Miura transformation (|6.ip . which to 
order A reads 

g x g- l = [VA] = \[JA]. (6.7) 
This imposes the following condition on g, 

{J,g x g~ 1 }=o, (6.8) 

which can be expressed as {S,g~ 1 g x } = 0, hence together with (|6.5p we have 

g~ l g x = -\s x S=:^ 1 \ (6.9) 

Inserting this in the expression for 5j 2 , and using the above identities, leads to 

S t2 =S XX - {S X S{I + S)) X = \[S,S XX \ = 2jf\ (6.10) 

which is the (generalized) Heisenberg magnet equation (see also Remark [5]). 

Remark 8. The ordinary Heisenberg magnet equation St = S x S xx is obtained from the 2x2 
matrix case by writing S = Ylk=i ^u^ki where ci, (72,(73 are the Pauli matrices, and setting 
*2 = -it- 



More generally, with the help of the Miura transformation (|6.ip the hierarchy equations 
resulting from (|6.3|) can be expressed solely in terms of S. (|6.ip implies 



1 {j.fe-.- w te-'}-«p*-*- w n-( a<,, T w) o) = (" 2A f w 

= -2AP[P,0][P,^- W P = -2AP 9l < ) - 1 ( 9l9 - 1 )_ w P, 



where we assumed that (I3.14p holds and used (|6.7p in the last step. Multiplication from the left 
by g~ l and from the right by g leads to 

\-\S,g-\g-g_ [X[ )} = -^X{I + S)S x g- 1 g XiAX[ g-_\ x] g(I + S), (6.11) 

where we used g~ x Vg = ±(I + S), QESJ, and S 2 = I. At order A 1 we obtain 

{5, g" 1 ^ - to)} = -5(1 + S)S x g- 1 g x (I + S) = i(J + S)S X 2 (I + S) = l -S x 2 {I + S), 
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hence, using (|6.6p and (|6.9p . 

{S, g~ l g t2 } = {S, g- l g xx ] + l -S x 2 (I + S) 



1 -{S X S) X + -AS x Sf \ + ^-S X \I + S) = ^-S x \ (6.12) 



With the help of ([6~l"0t ([63]) together with <^A2h implies 
1 3 

At order A 2 , flO) yields 

5*3 = Q«5t 2 - ^ + \{g" l g t2 + Or 1 ) xx g) (I + S)) . 

With the help of (|6,9p and (|6.13p . this can be arranged into the form 
St 3 = ^— — S XX S — ~^S X 2 — —S xx + — S XX (I + S) + — S X 2 (I + 5) 
Sxx + ^SJ = 27< 2 >. 
In a similar way, from (|6.3p and ()6.1ip we obtain 

9 1 9ta = —~^{SxxS) x — S X S XX — — S X 3 S =: J-^ 3 \ 

g 1 gt4 = ~^s x 2 s xx + — S X S XX S X + — s x s xxx s + — s xx s x 2 + — s xx 2 s 

_i_ ^ c c c i -9 4 U^ -■ tt(4) 
4 16 2 

and 

S tn =2F x n ~ l \ n = 2,3,4,5, 

which likely extends to all higher n G N. Of course, the expression for 5i n can be recovered by 
inserting the corresponding expression for g~ 1 gt n in (|6.5p . 

Remark 9. Conditions like (|6.8|f ^l originated from the use of the Miura transformation, and 
they are in fact needed to express the original hierarchy for the matrix variable g in terms of 5. 
The "mismatch" in the Miura transformation, leading to the restriction of the form of g, can be 
traced back to the fact that in the step from (|2.4p to (|2.3p we are dropping cohomological terms 
(see Remark 

Example 2. If m = 2, we write 

_i / k —a \ . I R a 

g = r I , hence g - 



ok/ kk — aa \ o k 



12 Via a transformation g H y ag, with a suitably chosen invertible block-diagonal matrix a, we can always 
achieve that (|6,8|) holds (note that the anticommutator of any matrix with J is block-diagonal) . Though this is 
a symmetry of the Heisenberg magnet equations (since it leaves S invariant), it is not a symmetry of the equations 
for g. 
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and thus 



S = g- 1 Jg 



1 f kk + aa 2ko 
kk — go V —2Ro —KK — aa 



The condition (|6.8p amounts to 
equations, we find that 

(kk — aa) x = 0. 



and kk x — ao x = 0. By adding these two 



(6.14) 



Remark 10. Using the Miura transformation (I2.3p . the linear equation (|3.2ip reads 

d^ = [(dg)g- 1 ]^ + 2zd^ 
(writing v = 2z), hence 

d^ = 2z[di> + g~ x {A.g)i>] 
in terms of 

$ = g~ 1 V. 

Evaluation of the linear equation for the bidifferential calculus given by (|3.7p leads to 

A- 1 ^ - A _1 ^_[ A ] (I - 2z\V) = 2z(gZ\ x {Pg)i>. 
Writing 

^ = ^ e -£»>i( 2 *p) B **, 

we obtain 

X -\^ - V>_ [A] ) = 2z(gZ 1 [x] Vg)^ = z(gZ\ x] g) (I + S)^. 

Expansion in powers of A yields 

ijj x = z(I + S)ip, 

4> t2 = [2z 2 {I + S)-zS x S]i;, 

^ 3 = [4z 3 (/ + S) - 2z 2 S x S + ^{2S XX + 3S X 2 S)] V, 



4z 3 (/ + S) - 2z 2 S x S + -{2S XX + 3S X 2 S)^ 
8z 4 (I + S) - 4z 3 S x S + z 2 (2S XX + 3S X 2 S) 



2 {^SxxxS ~ 2S XX S X — AS X S XX — 5S X S S) 



The first two equations constitute a Lax pair for the generalized Heisenberg magnet equation 
(with S 2 = I). 
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6.2 A class of solutions 

The following result is an analog to that in Section 13.31 (see also Remark 3 in pQ ) . It allows to 
generate solutions of (|2,4p from solutions of a linear system. 

Theorem 2. Let (f2,d,d) be a bidifferential calculus with $7 = A® f\(C N ) and A = Mat„, (£>), 
for some uq G N. For fixed n > no, let S G Mat(n, n, B) and A G Mat(m, m, B). Furthermore, 
let X G M&t{n,n,B) and W G Mat(m,re,Z?) satisfy the linear equations 

dX = (dX)S, dW = (dW)S, 



and also 



XS - RX = VZ, WS - AW = CX, (6.15) 

with matrices C,Z G Mat(m,n,B) ; -R G Mat(n,n,B) and V G Mat(n,m,B) ; satisfying 
dC = 0, dR = 0, dV = 0, dV = 0. 

Then 

g={WX- l V)-\ (6.16) 
provided the inverse exists, solves the {modified) Miura transformation equation (|2.3p . i.e. 

[d 5 - (d^A]^- 1 = d4>, (6.17) 
w'f/i some m x m matrix (f), and thus {by application of d) ako (|2.4p , z.el^l 

d{{dg-{dg)A)g- 1 ) = 0. (6.18) 
Proof. Using the Leibniz rule and the assumptions, we have 

d^- 1 = {dW)X~ 1 V - wx- 1 {dx)x- 1 v 

= (dW)SX~ 1 V - WX~ 1 {dX)SX- 1 V = (dW - WX- 1 dX)SX~ 1 V 

= d{WX~ l )XSX- l V = d{WSX~ 1 V) - W X~ l d{X S X^V) 

= d{Ag~ l ) - WX- l d[{RX + VZ)X- 1 V] = d(Ag~ l ) - g- 1 d(ZX' 1 V), 

and thus 

[dg - {dg)A]g- 1 = d{ZX' 1 V - gAg~ l ) . ■ 

Remark 11. The assumptions in Theorem [2] give rise to integr ability conditions. The latter 
are satisfied if 

dS = {dS)S, dA = (dA)A, dC = {dA)C, dZ = {dZ)S. 

In the following we exploit Theorem [2] for the bidifferential calculus given by (|3.7|) with some 
simplifications. We set A = and make the further assumption that S is d- and d-constant, and 
we write Z = UY, where U G Mat(m, n, B) is d- and d-constant and Y G Mat(n, n, B) solves 



13 This equation is invariant under right multiplication of g by any invertible d- and d-constant matrix that 
commutes with A. 
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dY = (dY")<S. This is motivated by the fact that then the first of conditions (|6.15p reduces 
to (13T23D . i.e. 

SK-KS = VU, (6.19) 

assuming that A d is invertible and using results from Section 13.31 in particular the defini- 
tions (|3.26p for K, U, V. Furthermore, we obtain 

W = W d + W S, S = e~^ J , £(S) = Y,S%, 

k>l 

and the second of conditions (16.15h yields C = W d SA^ 1 (which simply determines C) and, 
assuming that S is invertible, 

W = -W d SKS-\ 
The solution (|6.16|) of (|6.18|) (with A = 0) is then given by 

g- 1 = W d {l - SKS^Ejil - KSyW. 

Rewriting this as 

g- 1 =W d (l- SKS^S)^ + JSCS) (I - (KSfyW 

= W d (I - SKS'^S - SKS^SKS + KB) (I - (KH) 2 )~V, 

it can easily be decomposed into a part that commutes with J, 

{g' 1 ), = W d {l- SKS^EKS) {I - (KS) 2 )- V 

= W d V + W d {KS - SKjS^EKS^ - (K3) 2 )~V 
= W d V[l - US^SKS^I - (KSfyW] , 

and a part that anti-commutes with J, 

(g- 1 ) o = W d (KS - SK)S~ l E{l - (KS) 2 yW 
= -W d VUS~ l E(l - (KSfyW. 

Using our concrete form of J and «7, the matrices K, S, U, V have the form given in (13.32() . 
and we have 

This leads to 

»-' - ( —if 1 ) 

where 

K = (WV) [I + US^ZKZil - KEKE)~ l V] , 
R = (WV) [I + US^EKEil - KEKE)- l V] , 
a = -(WV)US- X E{I - KZKZy 1 ?, 
a = -{WV)US- X Z(I - KEKE)-^. 
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The only restrictions that have to be imposed on the matrices K, K, S, S, U, U, V, V result 
from (|6,19|) , They are 

SK + KS = VU, SK + KS = VU. (6.20) 

The solutions of the hierarchy for g obtained in this way also determine solutions of the generali- 
zed Heisenberg hierarchy. This is so because the solutions constructed above via Theorem [2] are 
actually solutions of the Miura transformation and our choice of matrix data via Proposition [T] 
ensures that (|3.14p holds (which we used in Section 16. ip . 



6.3 Reciprocal dual and combined dual AKNS hierarchies 

Elaborating the dual equation (|2.4p with the "reciprocal" bidifferential calculus determined 
by (|4.ip . instead of using that determined by (13. 7p . we simply obtain (|6.2p with g replaced 
by g" 1 . Again, we can combine the dual AKNS hierarchy and its reciprocal version, adopting the 
procedure in Section [5j New equations arise from the mixed parts, hence from evaluating (12.41) 
using the bidifferential calculus given by 

df=[VE x J]C + ^ l [E fl J}C, df = \- l [E x ,f]C + [VE ll ,f}C, 

which is a constituent of the calculus determined by (|5.ip . This results in 

[Pg-toPg- 1 - (g-viPg- 1 )-^] - m"^" 1 fo-MS" 1 - (W 1 )^] = °- 

To order X°fjP this is 

{gxg-^s+fatfPg- 1 ] =0, 

hence 

(gxg' 1 )^ = \ igg~ 1 - gg' 1 ) where g = J g J - ( 6 - 21 ) 

The Miura transformation between the combined hierarchies consists of a pair of Miura 
transformations, one for the original hierarchy and another one for the reciprocal. It results in 
the following two generating equations, 

*- 1 {g-g-[\ ] )g~ 1 ='P<f>-<t>-lx]'P, {Vg- g_ [x {P)g- 1 = X- 1 ^- 4>_ [x] ). 

In particular, this yields 

9xg- 1 = [V,<f>], <f>x = [P,g\g- 1 - (6.22) 

The formulas in Section T6.2I still generate solutions of the combined dual hierarchy and also of 
the Miura transformation (cf. (|6.17j) ). provided we extend the expression for £(S) used there to 



as) = Y,s k t k +Y, s ~ k tk- 

k>l k>l 



(6.23) 
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6.3.1 Sine-Gordon solutions 

If m = 2 and if g has the form 

/ cos(tf/2) -sin(0/2) \ _ M/2 ( -1 \ 

9 1 V sin(0/2) cos(tf/2) y /e ' \1 J 1 

where / is a function independent of x, then (|6.2ip becomes the sine-Gordon equation 

&xx = sin(i?). 

The function / drops out of equation (|6.2ip . As a consequence of the form of g, the condi- 
tion (|6.8p . which arose from the Miura transformation, is satisfied. (I6.22p requires the reduction 
conditions (|5.6p with e = — 1 and then reads 

q = -^x, fe = --sin(#), p s = -[1 - cos(tf)]. 

In order to generate solutions of the sine-Gordon equation (and more generally of the cor- 
responding hierarchy), we have to choose the matrix data in such a way that g has the above 
form. We set 

K = -K, S = S, U = -U, V = V, W = W. 

Then (16.20P reduces to a single Sylvester equation, SK + KS = VU. Setting t 2 k = ^2k = 0, 
k = 1, 2, . . ., (|6.23p has the property £(— S) = —£(£). As a consequence, we have 3 = e^ -5 ) = 
e — — ^ anc j thus 

K = r = a [l -US^EKE^ +(KE) 2 )~W], a = -a = -aUS^E^ + {KE) 2 )~ X V, 

where a = WV. We still have to ensure that k 2 + a 2 = f 2 , with some function / that does not 
depend on x. But since our procedure actually solves the Miura transformation (recall (I6.17P ). 
we already know that (|6.8p is satisfied, hence (|6.14p holds, which shows that k? + a 2 indeed does 
not depend on x. 

Example 3. Let n = 1, S = s £ M, and U = V = W = 1. Then the Sylvester equation 
SK + KS = VU is solved by K = Writing E = 2se~^ s \ where £(s) = Efe>o( s2fc+1 *2fc+i + 
s -2fc-1 t2fc+i)+£o with a constant we obtain k = tanh^(s) and a = sech^(s), so that k 2 +<t 2 =1. 
From a Ik = csch£(s) then follows the well-known kink solution $ = 2 arctan(csch£(s)) = 
4 arctan e^ s ). With n > 1 and real diagonal S we obtain multi-kink solutions. 

7 Conclusions 

We have shown in particular how a large family of solutions of matrix NLS equations, obtained 
in [2] with the help of general results of [I], extends to solutions of the corresponding hierarchies. 

Moreover, by a simple exchange of the roles of d and d, we obtained a "reciprocal" or 
"purely negative" counterpart of the AKNS hierarchy, which turned out to be the nonlinear 
part of the potential KP hierarchy. Combining the two hierarchies then gives rise to additional 
"mixed flows" . In this way we recovered in particular the short pulse equation and obtained an 
apparently new vector version of it (different from those considered in [421 144j ) . for which we 
presented soliton solutions in the 2-component case. 

Via specialization of the general Miura transformation to the bidifferential calculus studied in 
this work, we recovered a relation between the AKNS hierarchy and the "dual" hierarchy of the 
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generalized Heisenberg magnet model. As the first "mixed flow" of the dual hierarchy combined 
with its negative counterpart, with a certain reduction the sine-Gordon equation showed up. 

In this work we concentrated on a simple method, introduced in [I], to generate a class of so- 
lutions, parametrized by certain matrix data (essentially of arbitrary size) subject to a Sylvester 
equation. The largest part of the work in [2] concentrated on narrowing down a remaining 
redundancy in the matrix data that determine a matrix NLS solution. We expect that most of 
these results can be carried over to the cases treated in the present work. 
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